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ABSTRACT 

We report limits in the planetary-mass regime for companions around the nearest single 
white dwarf to the Sun, van Maanen's star (vMa2), from deep J-band imaging with Gem- 
ini North and Spitzer IRAC mid-IR photometry. We find no resolved common proper motion 
companions to vMa2 at separations from 3 — 45", at a limiting magnitude of J « 23. Assum- 
ing a total age for the system of 4.1 ± 1 Gyr, and utilising the latest evolutionary models for 
substellar objects, this limit is equivalent to companion masses > 7 ± 1 Afj„p (Toff ~ 300 K). 
Taking into account the likely orbital evolution of very low mass companions in the post- 
main sequence phase, these J-band observations effectively survey orbits around the white 
dwarf progenitor from 3 — 50 AU. There is no flux excess detected in any of the complimen- 
tary Spitzer IRAC mid-IR filters We fit a DZ white dwarf model atmosphere to the optical 
BVRI, 2MASS JHK and IRAC photometry. The best solution gives Toff = 6030 ± 240K, 
log g= 8.10 ± 0.04 and, hence, M = 0.633 ± O.O22M0. We then place a 3 cr upper limit of 
10 ± 2 Mj„p on the mass of any unresolved companion in the 4.5/im band. 
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1 INTRODUCTION 

Direct imaging of extra-solar planetary-mass companions to solar- 
type stars is complicated by problems of contrast and resolution. 
At the time of writing, no planet has been directly imaged around 
a solar-type star. An alternative solution to these difficulties is to 
target intrinsically faint stars instead, such as white dwarfs. Stel- 
lar evolution lends two huge advantages when searching for very 
faint companions: white dwarfs are up to ^ 10* times fainter than 
their main sequence progenitors, and the orbits of any planetary- 
mass companions that lie outside the stellar envelopes during the 
giant phases will expand outwards as mass is lost from the cen- 
tral star, increasing^ the projected separation by a maximum fac- 
tor Mms/A^wd yeans 1924). Thus, the problems of contrast and 
resolution are greatly reduced. The possible evolution of plane- 
tary syst ems in the po st-main sequence phase is discusse d in more 
detail by [ Duncan & L issauer ( 1998), Burleigh, Clarke & HodgkinI 



detail by 
( |2002|) ,Ie 



Debes & Sigurdssoij ( 120021) and lVillaver & Livid ( l2007t) 



The direct detection of such low mass companions to white 
dwarfs opens up the possibility for spectroscopic investigation of 
a previously unobserved class of object: evolved low-mass brown 
dwarfs and planetary-mass gas giants as low in temperature as 
left ~ 300 K. In contrast, the directly imaged planetary-mass 
5Afj„p) companion to the br own dwarf 2MASSWJ120 7 has 
the spectrum of a mid-L dwarf iChauvin et al. [ |2004l2005i) . be- 
cause it is still young (~ 10^ years old). The coolest known 
brown dwarf, ULAS J003402.77 - 005 206.7, has a temperature 
600K< Tefi < 700K jWarren et alj2007h and a spectral type T 8.5. 
The letter "Y" has been suggested for the next, cooler, spectral type 
dKirkpatrickl 1200 5). which might include evolved planetary-mass 
companions to white dwarfs. Alternatively, if no obvious spectral 
change triggers the use of a new letter, then the T classification will 
need to be extended beyond T 8.5. 

The idea of using white dwarfs t o find intrinsically 
faint, low mass companion s is not new. IProbstI il983h and 
iBecklin & ZuckermanI l ll988l) used the low luminosity of white 
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dwarfs to search for brown dwarf companions as near-infrared 
photometric excesses, and ind eed the latter achiev ed success with 
GD 165 B (L4). More recently. Iparihi et alJ Jioos') have conducted 
a comprehensive search for brown dwarf companions to several 
hundred white dwarfs but detected o n ly one new pair (G D 1400, 
L 6 - 7, [ Farihi & Christophe r! |2004 lOobbie et alJ l2005h . while 
iMaxted et al. (2006) and Bur leigh et al.l ( l2006h liave detected a 
close L8 companion to the white dwarf WD 0137— 349. No other 
detached substellar companions to white dwarfs are known. 

The most intriguing circumstantial evidence for the existence 
of old planetary systems that have survived to the final stage of 
stellar evolution comes from the discovery of metal-rich circum- 
stellar dust and gas disks aroimd a growing number of white dwarfs 



I Zuckerman & BeckliJl9"87l:lBecklin et alj2005l:lKilic et alJ200l 



2006l:lGansicke et alj2006l : lvon Hippel et al.l2007l : |jura et alj2007 



Kilic & Redfiel j2007 ). Spitzer mid-infrared spectroscopy has 
revealed that the dust disks are composed largely of silicates 
( iReach et aLll2005l : [jura et alj|2007l) . The favoured explanation for 
the origin of this material is the tidal disruption of an asteroid tha t 
has wandered into the Roche radius of the white dwarf jjurj2003h . 
perhaps through interaction with planets in a solar system that has 
become destabilized in the wake of the planetary nebula phase 
( iDebes & Sigurdssorj|20o3) . ljural J2006h suggests that at least 7% 
of white dwarfs possess asteroid belts. If that is indeed the case, 
it is likely that at least this number of white dwarfs also possess 
planetary systems. 

The recent discovery of a A/sini = 3.2Mi„p planet in a 1.7 AU 
orbit around the ex treme horizontal branch star V391 Pegasi by 
ISilvotti et al.l j2007h proves that such objects can survive red gi- 
ant branch evolution, strongly suggesting that it is simply a matter 
of ti me before a planet is d iscovered around a white dwarf. 

iBurleigh et al.l 1 I2OO2I) made predictions concerning the likely 
near-infrared brightness of putative resolved, giant planetary com- 
panions to nearby white dwarfs, based on their likely total ages and 
distances. In 2002 we initiated a programme to search for wide, 
spatially resolved very low mass common proper motion compan- 
ions to white dwarfs via direct imaging. In particular, we aim to 
find substellar companions with masses greater than a few A/j„p 
around white dwarfs within « 20 pc of the Sun. Such companions 
are expected to have near-IR magnitudes brighter than J ~ 23.5, 
commensurate with the expected sensitivity of an 8m telescope in 
a one hour exposure. 

We christened our project "DODO" - Degenerate Objects 
around Z)egenerate Objects. Preliminary results and progress re- 
ports have been published elsewhere ([C larke & Burleigh 20041 : 
iBurleigh. Hogan & Clarkel2006l : lHogan. Burleigh & Clarkel20^ . 
Here, we report our results for the nearest single white dwarf to 
the Sun in ou r sample, van Maanen 's star (vMa2, WD 0046+051, 
d= 4.41 PC. fPerrvman et alj |l997'). We combine two epochs of 
deep, ground-based J-band imaging from the 8m Gemini North 
telescope with Spitzer mid-infrared photometry to place limits on 
the masses and temperatures of any common proper motion and 
unresolved ultra-cool substellar and planetary-mass companions. 



1.1 vMa 2 and previous infrared observations 

vMa 2 was discovered serendipitously by Ivan MaanenI jl917l) in a 
survey for common proper motion companions to an unrelated star, 
HD4628. It has a cool, helium-rich atmosphere and strong reso- 
nance Ca II H & K lines in its optical spectrum, and is classified as a 
DZ white dwarf. The heavy elements are expected to sink below the 
photosphere on a timescale much shorter than the white dwarf cool- 



ing time. Their presence has been explained previo usly in terms 
of ep isodic accretion from the interstellar medium dDupuis et al.l 
1 19921) . but this scenario requires the accretion rate of hydro gen to be 
at lea s t two orders of mag n itude lower than the metals ( Wolff et al.l 
I2OO2I : iDupuis et alj 11993 : iDufour et ^l2007h . Alternatively, the 
presence of metals in DZs might be explained from accretion of 
cometary material or tidally disrupted asteroids and planets. How- 
ever, no DZ has been identified with an infrared excess due to dust 
emission, possibly because the diffusion timescale for heavy el- 
ements in cool white dwarf helium-rich atmospheres is long and 
the material may have been accreted as much as ~ 10® years ago. 
Nonetheless, DZ white dwarfs are a speculative candidate for host- 
ing old planetary systems and vMa2, as the nearest single white 

dwarf to the Sun, is an ideal target for su ch a search . 

Through analysis of Hipparcos data. lMakarovl j2004h claimed 
to have detected astrometrically a 0.06 ± O.O2M0 substellar com- 
panion to vMa2, with an orbital per iod of 1.57 years and a max- 
imum separation on the sky of 0.3". lFarihi. Becklin & MacintoshI 
(2004) carried out a search for this companion by direct imaging 
with adaptive optics in the mid-infrared L' band. They also looked 
for an unresolved companion as a near- and mid-infrared photo- 
metric excess in ground-based and ISO photometry. They refuted 
the existence of Makarov's companion, and of any excess emission 
due to dust, and placed a limit of TcS ^ 500 K on the temperature 
of any putative substellar companion. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Gemini NIRI ,7-band imaging 

Two deep J band images centred on vMa2 were a cquired using 
Gemini North and the Near-Infrared Imager (NIRI. iHodapp et al] 
I2OO3') with the f/6 camera (0.117" pixer\ field of view 120" x 
120") on 2004 December 25th and 2005 August 25th, in queue 
scheduled mode under programmes GN-2004B-Q-23 and GN- 
2005B-Q-19. The images were acquired with 60 s integrations at 
each position in a 54 point dither pattern, for total exposure times 
of 54 and 75 minutes respectively. The measured FWHM of stars 
in the images was 0.39" and 0.69" in the two epochs respec- 
tively. The data were reduced in the standard manner using the 
Gemini NIRI package in IRAF. A more detailed description of 
o ur data reduction method will be described in a later publication 
faogan et alfcoosi) . We note here that both these images were de- 
graded by intermittent 60 Hz electronic pick- up noise, which man - 
ifests itself as a diagonal herringbone pattern dHodapp et al.ll2003h . 



2.2 Spitzer IRAC mid-infrared observations 



Mid-infrared observations of vMa2 were o btained with the In- 
frared Array Camera (IRA CjFazio et alj2004b on the Spitzer Space 
Telscope (Werner et al. 2004j) on 2004 July 4th as part of Guaran- 
teed Time programme PID00033 (PI: G. Fazio). IRAC obtains im- 
ages in four channels with central wavelengths of 3.6/im, 4.5/im, 
5.8/im and 8.0/im. We downloaded the Spitzer IRAC data from 
the public archive as Basic Calibrated Data (BCD), which have 
been reduced and flux calibrated with the S14 version of the 
IRAC pipeline. The BCD images were corrected in two differ- 
ent ways as a check on our reduction method. Firstly, we cor- 
rected the individual images for array-location dependence, us- 
ing the correction images provided by the Spitzer Science Center 
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Table 1. Adopted parameters for van Maanen's star 
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" Hipparcos measurements, from our new fit to the optical, near-IR and IRAC mid-IR photometry (see section 4.2), estimated using evolutionary models 
appropriate for He-atmosphere white dwarfs with C/O core compositions , see http://ww w.astro. umontreal.ca~ bergeron/CollingModels/ derived from the 

initial-final mass relation of lOobbie et al.! i2006l) . " IWoojfigi^) 



(SSC), before carrying out the photometry on the individual cor- 
rected BCDs. Secondly, we corrected the BCDs for array-location- 
dependence using the coiTection images that had been divided by 
the pixel-distortion images (also provided by the SSC). We com- 
bined these corrected BCDs with dual-outlier rejection to produce 
a mosaicked image for each channel using the SSC mosaicking 
software, MOPEX dMakovoz & Marleaul llOOl) . By dividing the 
array-location-dependence images by the pixel distortion images, 
we avoid doubly correcting for pixel distortion when combining 
the BCDs with MOPEX. Both of the reduction methods gave the 
same results to well within uncertainties. 



3 DATA ANALYSIS 
3.1 J-band images 

Each Gemini NIRI J-band image was set to the World Co-ordinate 
System using stars in t he field with good 2MASS J-band detections 
dSkrutskie et"ai]|2006h . We used the SExtractor task to detect all 
point sources in the final stacked image with a flux ^ 3 cr above 
deviations in the local background, and computed the motions of 
these objects in RA and Dec between the two epochs (Figure 1). 
We also show in Figure 1 the 1 a and 3 u scatter of the distribution 
of all the measurements (excluding the white dwarf). 

The instrumental magnitude of each detected source was cor- 
rected to apparent magnitude _using stars in the field with 2MASS 
J-band fluxes dSkrutskie et alj|2006 ). To calculate the 3cr com- 
pleteness limit of these data, we inserted a total of 10,000 fake stars 
into each image at a variety of magnitudes from J = 19 — 24, and 
attempted to recover them using SExtractor. Figure 2 shows the 
fraction of implanted stars recovered as a function of their bright- 
ness at each epoch. The number of implanted stars recovered is al- 
ways less than 100% of those injected as some stars are lost behind 
or within the wings of other real or implanted stars. The fake stars 
were inserted at the same positions in both images. However, each 
implanted star may not necessarily be recovered in both images, 
especially towards the completeness limit, in which case it would 
not be available for a proper motion measurement. The curve with 
the star markers takes this into account. 



3.2 Spitzer IRAC mid-infrared photometry 

We performed aperture photometry on the mosaicked IRAC im- 
ages usin g standard IRAF tools. A full description is given in 
iBrinkworth et al. (2007). Briefly, we adopted an aperture size of 
2.5 pixels, and subsequently corrected to an aperture size of 10 pix- 
els, using corrections of 1 . 1 90, 1 . 1 99, 1.211 and 1 .295 for channels 
1 — 4 respectively. The IRAC fluxes of vMa 2 are given in Table 2. 
The errors consist of the quadratic addition of the standard devia- 



Figure 1. Proper motions of all point sources detected in both the NIRI 
J-band images, between 2004 December 25th and 2005 August 25th. The 
white dwarf's motion is clearly discriminated from all other objects in the 
field. The dashed circles represent the 1 a and 3 a scatter of the distribution 
of the proper motions, centred on the white dwarf for clarity. 
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Table 2. Predicted and observed Spitzer IRAC fluxes 





Predicted 


Observed 


Ao 
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(mJy) 


(mJy) 


3.56 


7.991 


8.030 ±0.262 


4.51 


5.370 


5.451 ±0.170 


5.76 


3.551 


3.611 ±0.126 


7.96 


2.038 


2.092 ±0.089 



tion of the fluxes in individual frames plus th e systematic error on 
the calibration of each filter as determined bv lReach et alj ( l2005h . 



4 RESULTS 

4.1 Limits on common proper motion companions 

Figure 1 shows the proper motions of all point sources detected 
in the NIRI images, between the two epochs. The white dwarf 
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Figure 2. The completeness limit of each NIRl observation is determined 
by injecting fake stars into each image and recording the fraction recovered 
by our source detection method against magnitude, as described in the text. 
We also show the fraction recovered in both images. 
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is clearly discriminated from everything else in the field, i.e. we 
do not detect any common proper motion companions to vMa2. 
We recover 90% of fake stars injected into both images at J ~ 
20.9 mags, and 50% at J ~ 22.7 mags (Figure 2). These sensi- 
tivity estimates can be translated into limits on the mass and tem- 
perature of any resolved companions, using appropriate evolution- 
ary models for subs tellar objects. We adopt the COND models of 
lBaraffeetal] ( l2003l) . which give the temperature and luminosity of 
substellar objects evolving in isolation, i.e. we assume no insola- 
tion from the white dwarf, and neglect any possible heating during 
previous evolutionary phases. The age of vMa 2 can be calculated 
by combining the white dwarf cooling age (3.17 ± 0.29 Gyr, see 
below) with an estimate of the main sequence l ifetime. Using the 
initial-final mass relation of iDobbie et al.l ( l2006h . we estimate that 
the progenitor of vMa2 (Mwd = 0.633 ± 0.022 Af©, see be- 
low) had a mass = 2.&Mq. We estimate the lifetime of such a 
main sequence star as ~ 0.9 Gyr, from the relationship of I Wood 
( 11992 '). Thus, we adopt a total age for vMa2 of 4.1 Gyr. At this 
age, and taking the Hipparcos-deriwed distance to vMa 2 of 4.41 pc, 
the 90% completeness limit (J ~ 20.9 mags) translates to a limit 
on the mass of resolved companions of ~ O.OO9M0 (~ 9i\fj„p, 
Tcff ^ 320 K), and the 50% limit (J « 22.7 mags) translates to a 
limiting mass of ^ O.OO7M0 (« 7Afj„p, T^s « 280 K). 

These mass and temperature estimates are of course model- 
dependent. For example, the white dwarf initial-final mass re- 
latio n is empirical and the subject of several current studie s 
(e.g. iDobbie et al]|2006l ; iKalirai et alj|2005l : iFerrario et aObOOSh . 
white dwarf evolutionary models are constantly being revised and 
improved, and the lifetimes of field main sequence stars are notori- 
ously difficult to estimate. At ages > 1 Gyr the evolutionary mod- 
els for substellar objects predict that they cool very slowly, and are 
relatively insensitive to age. Therefore, even if we adopt a very con- 
servative error of ±25% (±1 Gyr) on the total age of vMa2 then 



Figure 3. Mid infra-red spectral energy distribution of vMa2 from 3fim 
to 10/-tm. The observed Spitzer IRAC fluxes are shown by error bars. 
The DZ white dwarf model monochromatic flux is shown as a solid line, 
and the model fluxes averaged over the filter passbands are indicated by 
filled circles. Also sho wn for comparison are the 3.76/jm L' band flux 
from.Farihi et al.l j2004h and the ISO 6.75fim LW2 flux from lCharv et alj 
(1999 ). 
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the errors on the masses and temperatures of putative companions 
are small: 9 ± IMjup (300 ± 20 K) at the 90% completeness limit 
and 7 ± lAfj„p (also, 300 ± 20 K) at the 50% completeness limit. 

The minimum projected physical separation surveyed is given 
by the distance from the centre of its point spread function where 
the white dwarf is indistinguishable from the background. We esti- 
mate this as ~ 3", equivalent to ~ 13 AU at the distance of vMa2. 
The maximum projected physical separation surveyed is given by 
the field of view covered by both NIRI images. Although the field 
of view of the NIRI f/6 camera is nominally 120" x 120", the 
dither pattern used during the observations restricts the useable 
field of view to a maximum radius of 45", equivalent to ~ 200 AU. 
As mass is lost from the central star during the post main se- 
quence phase, the projected star-planet separation for planets out- 
side the red giant envelope will increase by a maxi mum factor 
Mms/Mwd = 2.6Mq/O.63M0 « 4 i |jeans|[l92i) . From this 
we estimate the range of orbital radii surveyed around the white 
dwarf progenitor as 3 — 50 AU. 

4.2 Limits on unresolved companions through Spitzer 
mid-infrared photometric excesses 

In order to place limits on the masses and temperatures of unre- 
solved companions to vMa 2, we compare the Spitzer IRAC pho- 
tometry to an appropriate white dwarf model atmosphere plus 
COND models calculated for the IRAC bandpasses (Figure 3, Ta- 
ble 2). Firstly, we fitted a grid of DZ white dwarf model atmo- 
spheres to the optical (Bessel) BVRI, near-IR 2MASS JHK and 
the IRAC fluxes. The best fit solution gives Tch = 6030 ± 240 K 
and, when combined with the parallax, log g = 8. 10 ±0.04. Within 
the formal u ncertainties, these va lues are consistent with those de- 
termined bv lDufouretai] ( |2007|) . which were based on fitting the 
same BVRI data but a different set of CIT JHK photometry, and 
did not include the IRAC photometry. The IRAC fluxes now give 
us a much better leverage on the atmospheric parameters. We then 
folded this model through the IRAC filter passbands to predict the 
IRAC fluxes (Table 2). Our new predicted model is consistent with 
the observed fluxes in all optical and infrared bandpasses, within 
uncertainties. From evolutionary models, the best-fit parameters 
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lead to estimates for the white dwarf mass of 0.633 ± O.O22M0 
and cooling age 3.17 ± 0.29 Gyr. 

By adding COND models to the predicted IRAC 4.5/im flux, 
again assuming an age for vMa2 of 4.1 ± 1.0 Gyr and adopting 
a distance of 4.41 pc, we find that in this band Spitzer could have 
detected a 10 ± 2iV/,„p, T^s f» 280 - 360 K planet at the 3a level. 



5 SUMMARY 

We have placed limits on ultra-cool substellar and planetary-mass 
objects around the nearest single white dwarf, vMa2, through a 
search for common proper motion companions and mid-infrared 
photometric excesses. The red giant pr ogenitor to vMa 2 had a max - 
imum radius of ~ lOOOii© (4.6 AU. [Hurley. Pols & Tbu^bOOOh . 
Therefore, taking into account the post main sequence evolution of 
the orbits of any companions, our J-band images have covered all 
orbits of substellar objects that originally lay beyond the maximum 
extent of the red giant envelope, up to 50 AU distance. The compli- 
mentary Spitzer IRAC photometry places slightly higher limits on 
the masses and temperatures of any spatially unresolved substellar 
companion whose orbit fortuitously lies along the line of sight to 
vMa2, or that currently has an orbit within ~ 13 AU of the white 
dwarf. These limits are significantly low er t han those re ported by 
previous studies, e.g. lFarihi et al.l ( l2004h and lFarihi et alj f2005). 

A variety of searches for planetary companions to white 
dwarfs are currently underway (e .g. .Burleigh et all 1 20061 : 
iDebes et al.ll200^ : iMullallv et allioOTh . F urther, extensive re sults 
from the DODO survey are in preparation jHogan et alj2008l) . 
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